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Kondo Effect of a Magnetic Ion Vibrating in a Harmonic Potential 

Satoshi YASHIKI*, Shunsuke KIRINO, Kazumasa HATTORI and Kazuo UEDA 
Institute for Solid State Physics, University of Tokyo, Kashiwa 277-8581, Japan 

To discuss Kondo effects of a magnetic ion vibrating in the sea of conduction electrons, a 
generalized Anderson model is derived. The model includes a new channel of hybridization asso- 
ciated with phonon emission or absorption. In the simplest case of the localized electron orbital 
with the s-wave symmetry, hybridization with p-waves becomes possible. Interesting interplay 
among the conventional s-wave Kondo effect and the p-wave one and the Yu- Anderson type 
Kondo effect is found and the ground state phase diagram is determined by using the numerical 
renormalization group method. Two different types of stable fixed points are identified and the 
two-channel Kondo fixed points are generically realized on the boundary. 

KEYWORDS: Kondo effect, impurity Anderson model, electron-phonon coupling, numerical renormaliza- 
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1. Introduction 

Recently, a new type of ionic structure, a network of 
cages filled or unfilled by guest ions, has attracted much 
interest in the research area of condensed matter physics. 
When the radius of filled ion is smaller than the size of 
cage, the guest ions vibrate with large amplitude but al- 
most without effect of interference with oscillations of the 
ions in the neighboring cages. It is expected that various 
unusual phenomena are possible due to such incoherent 
vibrations. Typical examples of this class of materials are 
/3-pyrochlore oxides AOs206(A = K, Rb or Cs), filled 
skutterudites RT4Xi2(R — rare earth or alkaline earth; 
T = Fe, Ru, Pt or Os; X = P, As, Ge or Sb) and type-I 
clathrate compounds. 

The family of /?-pyrochlore oxides, especially KOS2O6, 
is a good example where several qualitatively new phe- 
nomena are observed due to the anharmonicity of the 
cage potential^^: concave-downward temperature depen- 
dence of resistivity^-*, peak structure of NMR relaxation 
rate at potassium site at low temperatures'^-' , occurrence 
of superconductivity'*-' and the structural phase transi- 
tion which takes place below the superconducting tran- 
sition temperature.^) Concerning the resistivity and the 
NMR relaxation rate, Dahm and one of the authors 
considered a simple model for the anharmonic poten- 
tial which includes the fourth order term of the ion dis- 
placement. They have shown that the unusual temper- 
ature dependence of these quantities are explained by 
calculating the retarded phonon propagator within the 
self-consistent quasi- harmonic approximation.^' Detailed 
properties of the spectral function of the local ion are 
also investigated for various types of the cage potential 
including a double-well type.^' The more realistic model 
that the ion vibrates in the three-dimensional cage poten- 
tial with the third order anharmonic term specific to the 
tetrahedral symmetry of the /3-pyrochlore compounds is 
discussed by one of the authors and Tsunetsugu.^"^' It 
is shown that the third order anharmonic term has im- 
portant effect on the superconductivity and is also re- 
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sponsible to the first order structural transition below 
the superconducting transition temperature. 

Another interesting behavior is reported for the filled 
skutterudite compound, SmOs4Sbi2.^°"^''-' It shows a 
large specific heat coefficient 7, which is robust against 
magnetic fields. ^^-^ One possible scenario of this behav- 
ior is that a non-magnetic Kondo effect is realized by 
local vibrations of the guest ions.*^'^^-* Motivated by the 
unusual heavy fermion behaviors of SmOs4Sbi2, we have 
started a systematic study on the interplay between local 
electron correlation and electron-phonon coupling. 

Let us review the previous studies related to Kondo 
effect coupled with local vibrations. To include effects 
of the local oscillations, two different types of models 
have been considered, the breathing mode of the cage 
and the transverse modes of the cage, i.e., the relative 
displacement between the ion and the cage. 

Firstly, a theoretical model of a magnetic ion coupled 
with the breathing mode is known as Anderson-Holstein 
model. In this model the one-electron energy level of 
the impurity orbital is modified by the amplitude of the 
breathing mode. By applying Lang-Firsov transforma- 
tion^^-', it can be shown that the interaction induces an 
effectively attractive electron-electron interaction (for ex- 
ample see Ref. 20). If the electron-phonon coupling over- 
comes Coulomb repulsion, the electronic states of the 
impurity result in a quasi degenerate doublet, the empty 
and the doubly occupied states. If there is a matrix el- 
ement between the doublet, it is possible that the non- 
magnetic doublet works as a pseudo spin and induces 
Kondo effect. Actually, a recent study shows that mag- 
netically robust Kondo effect is possible even in relatively 
high magnetic fields when the effect of anharmonicity of 
the breathing mode is included. 

Secondly, a model of a non-magnetic ion which oscil- 
lates in the sea of conduction electrons was proposed 
by Yu and Anderson. In this model, the conduction 
electrons are scattered by the vibrations of the ion. For 
simplicity spins of conduction electrons are neglected. 
Yu and Anderson considered only the first order term 
of ionic displacement which produced the scattering pro- 
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cesses between the ,s-wavc conduction electrons to the 
p-wave ones. It is shown that the resultant effective po- 
tential for the ion displacement behaves like a double well 
potential when the electron-phonon coupling is strong. 
Therefore a new type of Kondo effect is expected because 
logarithmic divergence appears when there are doubly 
degenerate states connected each other through scatter- 
ing processes of conduction electrons. ^^"^^^ 

This paper is a comprehensive report on the effects 
of vibration of a magnetic ion coupled with spinful con- 
duction electrons. Some of the results have been pub- 
lished in the previous letter. ^''^ In §2 we construct a gen- 
eralized impurity Anderson model of the magnetic ion. 
When we ignore spin-orbit interaction, the generalized 
gradient theorem^^^ enables us to expand hybridization 
term with respect to the ion displacement. As a result 
phonon-assisted hybridization is naturally induced in ad- 
dition to the conventional one. In §3 and §4, we apply the 
numerical renormalization group (NRG) method^^- ■^^^ to 
the simplest case, where the impurity ion with an s-wave 
orbital vibrates in a one-dimensional harmonic poten- 
tial. From the analysis of energy spectra obtained from 
the NRG calculations, one can identify three different 
types of NRG flows at low temperatures which deflne 
the low energy fixed points. We obtain the phase dia- 
gram characterized by the fixed points in the parameter 
space of phonon-assisted hybridization and Coulomb in- 
teraction. The fixed point analysis is described in detail 
in the present paper. Various physical quantities which 
include the spin and charge susceptibilities and the mean 
square amplitude of the displacement are also discussed. 
Finally, conclusions and some future problems are sum- 
marized in §5. 

2. A generalized impurity Anderson model for a 
vibrating magnetic ion 

2.1 Expansion of hybridization term by ion displace- 
ment 

In this section we derive a generalized impurity An- 
derson model for a magnetic ion vibrating in a cage 
potential. When the ion is static, it has to be equiva- 
lent to the usual impurity Anderson model. It is as- 
sumed that the effects of vibrations are taken into ac- 
count through hybridization modified by the shift of the 
ionic position. Conduction electrons are assumed to be a 
non-interacting system and decomposed into the partial 
waves, 



fecr 



(2.1) 



fecr £m 



where s{k) is the dispersion relation and Cka (c^,^) is an 
annihilation (creation) operator of the conduction elec- 
trons, cemaik) {cl^„{k)) is defined from Cka (4^) by 
using the spherical harmonics Yem{k), 

Cemaik) = V^^Y;^{k)cka. (2.2) 



The dynamics of the ion with position vector Q is de- 
scribed by the following Hamiltonian, 



= -|^(Vq)^ + K. 



(2.3) 



where K;agc(Q) is the cage potential and M the mass of 
the ion. Hybridization modified by the ion displacement 
Q is expressed by the overlap integral between a plane 
wave (l/-\/r2)e''^ '^ and the localized orbitals of the mag- 
netic ion, where Q is the volume of the system. With- 
out the spin orbit interaction, £-wa.ve electron orbitals, 
''Pemi'i') = ue{r)Y(m{r), satisfy the following Schrodinger 
equation. 



2m 



+ Kucleus(r') - et I V^m(r) = 



(2.4) 



where the potential V;iucicus('') is assumed to be isotropic 
and m is the mass of electron. When the ion is shifted 
from the origin by Q, the relative coordinate from the 
ion to electron coordinate r is given by i — Q. Therefore 
the overlap integrals are expressed by 



V, 



{Q) = Ji= d^re-^'-^iJemir - Q), (2.5) 
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where <? is a coupling constant. With these integrals, hy- 
bridization terms modified by the ion displacement are 
given as 

^hyb = J2J2 [^(fe;^-) (QHjema + h-c] , (2.6) 

ka m 

where fema {flma) ^'^^ annihilation (creation) op- 
erators for the localized £-wa,ve orbitals. We expand the 
overlap integrals with respect to the ion displacement Q. 
In this paper we focus on the first order terms of Q. The 
hybridization terms expanded within 0{Q) are given as 
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E 



ik)f, 



ma 



+ h.c. 



(2.7) 



where V^, V^'j!^'™ ^ and qi,,„ are defined in Appendix. It is 
obvious that the rotational symmetry 5*0(3) is preserved 
in the generalized hybridization term (2.7). 

Under the simplest assumption that the ion has an s- 
wave electron orbital, the hybridization terms (2.7) are 
rewritten based on the Cartesian coordinate system as 



hyb 



y2{voci{k)u+h.c.} 



ka 



+ E E {VicliAk)f.Qi+^-<^-}, (2.8) 

ka i=x,y,z 

where coo- (ci,^^) are the annihilation operators for the s 
(pj)-wave components of the conduction electrons. It is 
verified easily that 1//°'*^' = = Vi because of the 

rotational invariance. 
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For simplicity, wc study the situation that the ion vi- 
brates along one direction, say x-direction. Using a one- 
dimensional harmonic potential (l/2)Ma;^(3^, we get the 
Hamiltonian for this study, 

H =Hc + Hhyh + fflocah (2-9) 

H,=J2<k){cl{k)coAk)+clikMk)}, (2.10) 

Hhyb=^{l^OcL(fc)/a + ll.C. 

+ Vicl{k),U{a + a^) + h.c.}, (2.11) 
JJiocai =f>oja^a + efY,PaU + Uflftflh, (2.12) 

cr 

where Vi/y/2 is set to Vi and a (a^) is the annihila- 
tion (creation) operator of the phonon. For the localized 
electron orbital, the annihilation (creation) operator is 
expressed by fa (fa) and e/ is its energy given by eq. 
(2.4) and U the Coulomb interaction. When we restrict 
the isotropic degrees of freedom of the ionic vibrations 
to one particular direction, the SO{3) symmetry is low- 
ered to the discrete inversion symmetry; r ^ —r and 
Q^-Q. 

2.2 Comments on the symmetries 

There are three relevant quantum numbers for the 
present model (2.9) — (2.12). It is obvious that the to- 
tal electron number A^*°' and the z component of total 
spin Sl°^ are conserved. The last one is related to the 
inversion symmetry mentioned in the previous subsec- 
tion. The operation of the inversion V is described in the 
second-quantized notations. 



VfaV-' = fa, 

rcaAk)V-' = coa{k), 

Vcia{k)V-' = -Cia{k), 

VqV-^ = -q. 



(2.13) 
(2.14) 
(2.15) 
(2.16) 



We characterize the inversion symmetry by the total par- 
ity P which is defined by the sum between the number 
of p-channel conduction electrons Np and the number of 
harmonic phonons A^ph, namely P = Np + iVph = or 1 
(mod 2).^°' This is based on the fact that, within 0{Q), 
phonon-assisted hopping between the impurity orbital 
and the p-channel conduction electrons has to accom- 
pany absorption or emission of one phonon. This quan- 
tum number plays a crucial role in the analyses of the 
energy spectra. 

In the original model used by Yu and Anderson (YA 
modcl)^^\ it is assumed that a non-magnetic ion, which 
is trapped in a one- dimensional harmonic potential, vi- 
brates in spinless conduction electrons. The first order 
process of the displacement induces scattering processes 
between the s-wave and p-wave conduction electrons. No 
localized electron orbital is introduced in the model. The 



original YA model is written as 

^YA = J2{^o{k)cl{k)co{k) + ei{k)c\{k)ci{k)} 



+ y,{yAik)ci{k'){a + a^) + h.c.} 



k,k' 



+ huia^a, 



(2.17) 



where co(fc) (ci(/s)) are the annihilation operators for the 
spinless s-wave (p-wave) conduction electrons with the 
dispersion relation eo(fc) (ei{k)). 

In general, eo{k) and ei{k) are different. However if we 
assume that eo(fc) and ei(fc) are the same, an artificial 
symmetry comes into existence. The symmetry, which 
we call channel symmetry, is the invariance under the 
following transformation I, 



Ico(fc)I-^ =Ci(fc), 
Ici(fc)I-^ =co(A;). 



(2.18) 
(2.19) 



With the channel symmetry, the Hamiltonian (2.17) can 
be decomposed into two detached parts by introducing 
the bonding and antibonding orbitals, 

co(fc) + ci{k) 



Ce(fc) 



Co(fc) = 



V2 ' 

Co(fc) - Ci(fc) 

^/2 ' 



(2.20) 



(2.21) 



and there is no matrix element between the two potential 
minima in contrast to the implicit assumption made by 
Yu and Anderson. Therefore the breaking of the channel 
symmetry is needed for the realization of Yu and An- 
derson type Kondo (YAK) effect. Matsuura and Miyake 
showed that potential scattering processes of the second 
order of the ion displacement play such a role."^^^ 

These discussions are also valid for the spinful fermion 
case. As for the present Hamiltonian (2.9), the introduc- 
tion of the localized impurity orbital serves as a source of 
breaking of the channel symmetry. Therefore the Hamil- 
tonian (2.9) at U = is one realization of the YA model 
extended to the spinful case. 

3. Numerical renormalization group approach 

3.1 Mapping to Wilson chains 

It is well known that for the ordinary impmity Ander- 
son model the numerical renormalization group (NRG) 
method is a powerful one to calculate energy spectra and 
various physical quantities with high accuracy.^''' We 
apply the NRG method to the present model by dis- 
cretizing the two continuous conduction bands in the 
logarithmic energy scales characterized by A. Both the 
,s-wave and p-wavc conduction bands can be discretized 
in the same way because the dispersion relation £(fe) is 
assumed to have no angular dependence on fc. The dis- 
cretized Hamiltonian is given by the two one-dimensional 
half-chains corresponding to the s-wave and p-wave con- 
duction bands, which may be called Wilson chains. 
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Fig. 1. (Color online) (a) Graphical representation of the Hamil- 
tonian (3.1). The wavy line represents phonon-assisted hop- 
ping while the straight lines usual hoppings. (b) The YA model 
mapped to the Wilson chains. 



N-1 
n—O.a 

- E [^4../. + y^v\ja{a + at) + h.c] 



^(1) V ftf -1 




(3.1) 



Here, for a constant density of states p = 1/2D with the 
band width of 2D the 7i-th hopping matrix element ^„ is 
given by 



l-A- 



2n-3 



(3.2) 



All the parameters with tilde are multiplied by the con- 
stant factor 2/{D{l + A^^)}. From now on we treat the 
symmetric case, 2e f + U = 0, for simplicity. The Hamil- 
tonian (3.1) is block-diagonalized and each block is char- 
acterized by the set of quantum numbers, N^°^, 5'°* and 
the parity P as discussed above. We set various param- 
eters for NRG calculations as follows; A — 3.0, cutoff 
number of phonon excitations 50, band width D = 1.0 
and M = 15000 states are kept at each NRG step. 

A graphical representation of the Hamiltonian (3.1) 
is shown in Fig. 1(a). For comparison, the YA model 
mapped to the Wilson chains is also shown in Fig. 1(b). 
From these pictures, we understand that within the first 
order of ion displacement there is the channel symme- 
try in the original YA model unlike the present model 
when the same hopping matrix elements are used for 
both channels. 

3.2 Two-site problem 

To investigate the effect of multi-phonon excitations, 
we consider the two-site Hamiltonian consisting of / and 
Po sites in the Hamiltonian (3.1), 



Ho. 



(1) 



(3.3) 



-Vi EK./- + h.c.}(a + at) 

+ hwa^a, (3.4) 



(3.5) 



For U — the two-site Hamiltonian is solved exactly 
by using the canonical transformation called Lang-Firsov 
transformation,"'^^'' e''^ with 



S = iAE(-4^a + BlB,){a - at). 



(3.6) 



where A(j and are the annihilation operator for the 
antibonding and bonding orbitals, A„ = (1/\/2)(/(t — 
Po,a) and B„ = (l/\/2)(/o- + po.cr). The coupling con- 
stant is given by A = —Vi/hio. Then the transformed 
Hamiltonian is written as 



-iS^(O) iS 



= -^|E(-4^. + 5i5^)| +nwaU. (3.7) 

We obtain the doubly degenerate ground states |L) and 
|R), which may be called polaron doublet (PD), with the 
ground-state energy i?pD, 

\-l)=b\bI\^^c)ji). 



with \t) 
|R)=44|vac)Jr)^^ 

with k)p,= e^^(— ^)|0)p. 



PD 



4lf 



(3.8) 



(3.9) 
(3.10) 



These two degenerate states are characterized by the ex- 
pectation values of the ionic position q, 

{L\q\L)= -{R\q\R)=2V2\, (3.11) 

16A2 + 1 



(L|g^|L)=(R|g2|R) = 



(3.12) 



(3.13) 
(3.14) 



One can see that the PD is formed both for spinful and 
spinless cases and it plays a similar role as the localized 
spin for the usual spin Kondo effect. 

In the following discussion, it is important to note that 
the electron parts of the doublet are written as 

4^Il-a^)cl=|LS)e,+ |PS>eP 

44|vac)^j=|LS)^,-|PS)^j, 
where 

|LS),.:^:^Ii|M^|vac), (3.15) 
represents a local singlet with even parity while 

|PS),^ ^^'^°-%^'^"-^ |vac), (3.16) 

a pair singlet with odd parity. Similarly, the phonon parts 
of the doublet are written as 



l^)ph = leven) , +|odd) 



'ph' 



(3.17) 
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Fig. 2. (Color online) Schematic one-particle energy levels of a 
single Wilson chain with (a) an odd number of sites and (b) an 
even number of sites. 



(3.18) 



l'^>ph = |even)ph"|odd)pj^, 

where |even)pj^ (|odd)pjJ is a state which is a hnear com- 
bination of states with even (odd) number of phonons. 
The parity of |even)pj^ (|odd)pj^) is even (odd). Clearly 

|L) and |R) are not eigenstates of the total parity, P. 

The degeneracy of the PD is lifted by Coulomb inter- 
action -ff2-sLc' which may be treated by the perturbation 
theory. In fact energy shift from i?pD is obtained by the 
following eigenvalue equation. 



det 



/Ai?PD-(L|4ielL) 



V -(R|4ie|L) 



-(L|4ie|R> 



=det 



A£:, 



PD 



U -8X^ 
4 ^ 



U 
4 



AEpo - X , 



= 0. 



(3.19) 



Then within the first order perturbation of U, the eigen- 
states and eigenvalues are written as 



l±) 



|L)±|R) 



V2 



U 



E^^ = EpYi + -\l± e-«^^ 



(3.20) 
(3.21) 



It shows that Coulomb interaction works on the PD like a 
transversal magnetic field for a localized spin. The energy 
gap between | -I- ) and | — ), Apd = (C//2)e~*'^ , is small 
when the electron-phonon couping A is strong, of the 
order of unity. Unlike |L) and |R), | -f ) and | — ) are the 
eigenstates of the total parity. The parity of | + ) is even 
{P = 0) and that of | - ) is odd {P = 1). 

4. Results of the NRG calculations 

4-1 Classification of fixed points from energy spectra 

To begin with, we summarize the NRG calculations for 
the usual impurity Anderson model. After the map- 
ping to a Wilson chain, the discretized Hamiltonian is 
composed of the impurity site connected to the single 
semi-infinite Wilson chain, which we will write as / site 
and Si {i = ^ n) sites. The number of the NRG steps n 
corresponds to the energy scale of the system. When n is 
sufficiently large, energy flow approaches to the unique 
low-temperature fixed point. 

The stable fixed point corresponds to the unitary limit 
of Kondo effect. In the limit, the low-energy spectrum is 
isomorphic to that of the Hamiltonian consisting of two 



Fig. 3. (Color online) Graphical representations of the two types 
of the energy spectra. 



parts, the singlet bond between / and so sites and the 
remaining free Wilson chain from si to Sn- 

Concerning analysis of the fixed point, it is essential 
to note that the low-energy spectrum depends on even 
or odd of n. Eigenstates of the Wilson chain are speci- 
fied by the set of quantum numbers, [7V*°*, 5'*°*]. For a 
tight binding Hamiltonian with only off-diagonal matrix 
elements, single-particle energies are symmetric with re- 
spect to positive and negative energies. When n is odd, 
the middle one of eigenenergies must be zero. At half- 
filling Fermi energy should be also zero. Therefore, the 
ground states have four-fold degeneracy which are spec- 
ified by the quantum numbers, [rt — 1, 0], [n, ±i] and 
[n -|- 1, 0]. For an even n, there is no eigenvalue at E — 
and therefore the ground state is unique with the quan- 
tum number [n, 0], see Fig. 2. 

Now let us discuss the fixed points for the present 
model. After n steps in the NRG, the system consists 
of 2n -|- 3 sites, the localized / orbital and the two Wil- 
son chains, sq, si, • • • , s„ and po, pi, • • • , Pn- Our nu- 
merical results by the NRG method show that there are 
only two stable fixed points: s- and p-type. The s-type 
fixed point schematically shown in Fig. 3(a) is charac- 
terized by the energy spectrum which is isomorphic to 
that of the system composed of tightly bound f — sq 
singlet and two independent Wilson chains whose length 
is different, ri sites for the s-chain and ri + 1 sites for 
the p-chain. It is important to include the parity as a 
quantum number, [iV'°*, 5''°*, P]. Concerning the s-type 
fixed point the phonon part does not play any impor- 
tant role and the total parity P is simply determined by 
the number of electrons in the p-channel. From Fig. 3(a) 
it is straightforward to see that the quantum numbers 
of the degenerate ground states are those given in Ta- 
ble I depending on even or odd n. As an example, we 
show in Table II a part of the low energy spectra of the 
s-type fixed point which correspond to the one particle 
excited states. The other fixed point defined as p-type is 
depicted in Fig. 3(b). The energy spectrum at low en- 
ergies is isomorphic to that for the system with isolated 
f — Po singlet bond and two independent Wilson chains, 
n + 1 sites for the s-chain and n sites for the p-chain. One 
important non-trivial fact obtained by the present NRG 
calculations is that the quantum number for the part of 
f — Po singlet is always given by [2, 0, 1], odd parity in 
the same way as the two-site problem discussed in §3.2. 
Therefore the quantum numbers of the ground states for 
the entire system are given by those shown in Table I. 
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NRG step n 


even 


odd 


fixed point 


s-typc 


[2n + 2, 0, 0] 
[2n + 3, ±|, 1] 
[2n + 4, 0, 0] 


[2n + 2, 0, 0] 
[2ri + 3, ±|, 0] 
[2n + 4, 0, 0] 


p-typc 


[2n + 2, 0, 1] 

[2n + 3, ±|, 1] 
[2?i + 4, 0, 1] 


[2n + 2, 0, 1] 
[2n + 3, ±|, 0] 
[2n + 4, 0, 1] 



Table I. Sets of [Af'°S S*°S P] of the ground states at the n-th 
NRG step depending on the types of the fixed points. 



s-type fixed point 





[1, 0, 




[1, 0, 1 




0th 





1 


0.800048 


2 


1st 


1.600095 


3 


2.400143 


2 


2nd 


1.695754 


2 


2.495802 


6 


3rd 


3.200191 


1 


2.997488 


2 



Table II. The low energy spectra of the s-type fixed point at an 
even number of the NRG steps for the sets of the quantum num- 
bers, [1, 0, 0] and [1, 0, 1]. In this table, Af*°* is counted from the 
total number of sites. The energy eigenvalues and their degener- 
acy factors are tabulated up to the third excited states. The low 
energy spectra of the p-type fixed point can be obtained by the 
inversion of the total parity P. 

Lastly, we discuss the unstable fixed points which sep- 
arate the regions of the s- and p-type fixed points. Here, 
we briefly review the 2-channel Kondo model. The model 
consists of a local moment coupled antiferromagnetically 
with two types of conduction electrons. The NRG 
calculations show that 2-channel Kondo effect is realized 
when the two antiferromagnetic coupling constants, Ja 
and Jb, are perfectly equal. '^'^^ The energy spectra ob- 
tained from the NRG calculations are summarized in 
Refs. 34 and 35. The low-energy spectrum of the 2- 
channel Kondo model with Ja — Jb shows S0{5) rota- 
tional symmetry, in the sense that the degeneracy factors 
of the eigenstates are subject to dimensions of its irre- 
ducible representations. 

Returning back to the present Hamiltonian, it is con- 
firmed that the spectra of low-energy eigenvalues on the 
unstable fixed points are identical to those obtained for 
the 2-channel Kondo model with Ja — Jb, see Table III. 
Therefore, we describe these unstable fixed points as the 
2-channel Kondo (2-chK) ones. 

4.2 Phase diagram of the generalized Anderson model 
for a vibrating magnetic ion 
From the energy spectra obtained by the NRG cal- 
culations, we determine the ground-state phase diagram 
which is shown in Fig. 4. The parameter space is divided 
into two distinct areas. The left part is characterized by 
the s-type energy spectrum and the right part by the p- 
type one. On the boundary, the energy spectrum is of the 
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Table III. Correspondence relation of the low energy eigenstates 
of the 2-chK fixed point between the present model and the 2- 
channel Kondo model at an even number of the NRG steps. In 
this table, Q is defined by the difference between N'^°^ and the 
total number of sites. Spin degeneracy is not counted in this 
table. The right part of the table is from Table 12 in Ref. 35 
with spin S and channel spin S^h for the 2-channel Kondo model. 
*There are doubly degenerate states. 




0.14 0.16 0.18 0.20 

Vi 

Fig. 4. (Color online) The phase diagram in the parameter space 
of Vi and U for Vq = cu = 0.2. On the boundary, the energy 
spectrum is of the 2-chK type. The left region is characterized 
by the s-type energy spectra and the right one by the p-type 
one. In the region between the green and blue dashed lines, the 
plateau of Si^np at log 2 disappears. 

2-chK type. On the line oiU — 0, the s-type energy spec- 
trum always appears. Therefore the line of 2-chK fixed 
points does not merge together to the U — line but 
approaches it asymptotically. 

The NRG calculations reveal that physical properties 
of the system at finite temperatures are sometimes quite 
different even if two different sets of parameters belong 
to the region of the same fixed point at zero temperature. 
Therefore, in addition to the analysis of the low temper- 
ature energy spectra, we have to investigate finite tem- 
perature behaviors of physical quantities at each point of 
the parameter space. 

First, we focus on the left side of the 2-chK line. On 
the line of U — this model is equivalent to the spin- 
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Fig. 5. (Color online) The line of level crossing in the parame- 
ter space of Vi and U obtained from the three-site two-electron 
problem for Vq = u) = 0.2. The color grade corresponds to the 
magnitude of log | — Eq \ , where Ea and Eo are the ground 
state energies with even and odd parity, respectively. In the left 
part of the deep blue line, Ee is lower than Eo- 



ful YA model. In the weak electron-phonon coupling 
regime main source of kinetic energy comes from the hy- 
bridization with the s-wave conduction electrons and the 
phonon-assisted hybridization with the p-wave channel is 
a weak perturbation. In this paper we characterize this 
regime by the renormalized Fermi chain (RFC). When 
Vi increases, the PD is formed and the degeneracy is 
lifted at low temperatures by the Yu and Anderson type 
Kondo effect, expressed as s-YAK in Fig. 4. Now, we 
gradually increase Coulomb interaction U , starting from 
the parameter region where temperature dependence of 
physical quantities is of the RFC type. The impurity be- 
gins to behave as a localized spin and then the spin is 
screened mainly by the s-wave conduction electrons. This 
regime is called s-channel Kondo (s-chK) in Fig. 4. In the 
parameter region between RFC and s-chK, the essence of 
the physics can be understood based on the local Fermi 
hquid theory. 

Next, we discuss the right part of the phase diagram. 
Concerning the strong U regime, it is easy to understand 
that the local moment is screened mainly by the p-wave 
conduction electrons, which is called p-channel Kondo 
(p-chK). When U is not important compared with Vi, 
the picture of the PD is valid, while the temperature 
dependence of physical quantities is different from that in 
the s-YAK parameter region in addition to the difference 
in the energy spectra. We denote this region as p-YAK in 
Fig. 4. In the following subsections, we will discuss details 
of physical properties for various parameters, including 
the difference between the s-YAK and p-YAK. 

4-. 3 Three- site problem 

At this point, it is instructive to consider the three- 
site problem, which consists of /, so and po sites. The 
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Fig. 6. (Color online) U dependence of the difference between E^ 
and Eo for the three-site problem with two electrons: Vq = u) = 
0.2 and Vi = 0.174. 




Fig. 7. (Color online) (a) U dependence of litc-cP and |^o-cP 
obtained for the three-site two-electron problem with Vo = u) = 
0.2 and Vi = 0.174. (b) U dependence of (g^)^ and (g^)^, which 
are the expectation values of the square of the ion displacement 
calculated from ^ovcn and 'I'odd- 
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Hamiltonian is given by 

-site — 

+ Vi'^{plfa + h..c.}{a + a))+njjja)a 

+ ^(E 4/^-1)'' (4-1) 

where the site sufBx " 0" is suppressed. It is easy to obtain 
the exact eigenvalues and eigenstates by using the exact 
diagonalization method. In the following examples the 
number of phonon excitations is kept up to 50 and Vq 
and iij are fixed to 0.2. 

Figure 5 shows the magnitude of log \E^ — Eo\ in the 
two-dimensional parameter space of Vx and i7, where E,, 
{Eo) is the ground state energies with even and odd par- 
ity. In the figure there is a line of level crossing between 
Eg and Eo- In the left (right) part of the line;, E^ (Eo) 
is lower than Eq {Eg). Figure 6 shows U dependence of 
the difference between Eg and Eo with Vi fixed to 0.174. 
Clearly there is similarity between Fig. 4 and Fig. 5. 

From the view point of parity, every component of 
wave functions is classified into four groups, VPe-e, ^e-o, 
^o-e and ^"0-0- for example, means a linear combi- 

nation of the components whose electron parts have even 
parity and phonon parts odd parity. Clearly an eigenstate 
with even total parity is written as 

^even ^e-e ~t~ ^o-O) (^•^) 

while one with odd total parity as 

*odd = *o-e + *e-o. (4.3) 

Concerning the three-site two-electron problem, the 
bases for a given parity are listed in Table IV. 

Figure 7(a) shows U dependence of the amplitudes, 
l^'c-oP and |*o-eP, for Vi = 0.174. Because the total 
wave functions are normalized, |^'o-oP = 1 — l*e-eP and 
|*e-oP = 1 - l^'o-eP- In Fig. 7(b) we show the expecta- 
tion values of the square of the ion displacement, (q^)^ 
and (q"^)^, which are calculated from \I/even and ^'odd) 
respectively. 

We start from the two-site problem for the non- 
interacting case, namely Vq = U ~ 0, where the ground 
states |-|-) and | — ) or |L) and |R) are doubly degenerate. 
When Vq is introduced, we expect that the kinetic energy 
favors the configurations which belong to 4'c-o- Actually, 
the amplitude l^'o-eP shown in Fig. 7(a) is larger than 
0.5 at U = 0. As we have discussed in §3.2 the Coulomb 
interaction U favors the statc^ with odd total parity. The 
competition leads to the level crossing at small U. 

In the narrow region of Vi , the ground state makes the 
second level crossing at larger U. In the strongly corre- 
lated regime, a local moment is formed in the /-orbital, 
which should be eventually screened by the conduction 
electrons. The expectation values of square of the ion dis- 
placement are suppressed and become close to the ampli- 
tude of zero point fluctuations for both Steven and Vfodd as 
shown in Fig. 7(b). In this situation, the phonon parts of 
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1 


electron 


/t5||vac)^j, 4/il™c)^i 


6-|p||vac)^j, p|s||vac)^j 


phonon 


■"),,], (" = I'vcn) 


■"-)p]i (r/ = odd) 



Table IV. Bases for electron and phonon parts of the three-site 
two-electron problem classified by the parity P. 

the wave functions are dominated by components with 
even phonon numbers. This fact is clearly seen in Fig. 
7(a). Therefore the total parity is predominantly deter- 
mined by the nature of singlet formation, even parity 
for the f — s singlet while odd parity for the f — p sin- 
glet. The upper critical U shown in Fig. 6 corresponds to 
the point where the two screening mechanisms compete 
in the same way as the 2-chK fixed point shown in the 
phase diagram for the bulk system. Fig. 4. 

The line of level crossing shown in Fig. 5 is qualita- 
tively very similar to the phase boundary obtained by 
the NRG calculations, Fig. 4. We would like to point out 
that the three-site Hamiltonian corresponds to the initial 
step of the NRG procedure and the ground state of the 
two-electron problem has the different parity depending 
on the nature of the singlet bond as discussed in §4.1, 
see also Table I. 

4-4 Strong correlation regime {U ^ w ^ T) 

Firstly, we discuss physical properties in the strong U 
regime where the relation, U ^ w ^ T, holds. In this 
regime multi phonon excitation processes are not impor- 
tant. Because of the strong U, the localized orbital be- 
haves as a local moment which is coupled with the s-wave 
and p-wave conduction electrons. With lowering temper- 
ature, the local moment is screened by the conduction 
electrons and the main screening channel changes de- 
pending on relative strength between Vq and Vi. In the 
same way as the conventional impurity Anderson model 
is mapped to the spin 1/2 Kondo model, we consider 
the projected Hilbert space where the vacuum and dou- 
bly occupied states of the impurity orbital are prohibited. 
We treat both of the conventional hybridization term Vq 
and the phonon-assisted one Vi in the Hamiltonian (2.9) 
as perturbation terms. After some calculations, a gener- 
alized Kondo model is obtained as 

i?Kondo=E e{k)cl^{k)ce„{k) + fuva^a + Hs.d, 

fecr £=0,1 

= - EE{-^O.o4a(fc)M„^CoMfc') • S 
fefe' a/3 

+ -^1,0 9[cL(fc) M^^co^fc') • S + h.c] 

+ Ji,i q^cUk)[a] c,p{k')-s], (4.4) 
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Fig. 8. (Color online) Temperature dependence of entropy in the 
strong correlation regime. The used parameters are Vq = <^ = 0.2 
with the symmetric condition U = —'^£f = 2.0. (a) For Vi < 
0.1742, the s-chK effect is observed, (b) When Vi is larger than 
the critical value the local moment is mainly screened by the 
p-wave conduction electrons, corresponding to the p-chK effect. 
At Vi = 0.1742, we observe signature of the 2-channel Kondo 
effect. 



where S — /^[cj^^/^ with tr being the vector of Pauh 
matrices and q ^ {a + a^)/\/2. The coupUng constants 
Jo.o, Ji,o aiid Ji,i are written as 



^^0,0 =Vq s 



1 



1 



Ji^o =Vo{V2Vi) 



1 



U 



u 



Ji,i =2V^ 



1 



u 



1 

.£/ 

The term proportional to Ji_o corresponds to scatter- 
ing processes from the s-wave conduction electrons to 
the p-wave ones and vice versa. Similarly, the term pro- 
portional to Ji_i describes scattering processes within 
the p-wave conduction electrons. In the strong correla- 
tion regime, effects of electron-phonon coupling are sup- 
pressed by Coulomb interaction. Therefore we may re- 
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Fig. 9. (Color online) Temperature dependence of spin suscepti- 
bility Txa in the strong correlation regime, = = 0.2 with 
the symmetric condition V = —2ef = 2.0. When Vi is relatively 
small, Txs increases close to 0.25 which corresponds to the local 
moment limit. With increasing Vi , phonon-assisted hybridization 
suppresses the enhancement of Txs- 



gard the effects of phonon emission or absorption pro- 
cesses as a renormalization of the coupling constant by 
taking average over the phonon degree of freedom, 

kk' a/3 



(4.5) 



where used is the fact that (g), the thermal averages of q, 
vanishes. Even if the electron-phonon coupling is weak, 
(g^) is not zero due to the zero point fluctuations. 

In Figs. 8 and 9 we use a set of parameters, Vq ~ 
u! — 0.2 with the symmetric condition U = —2ef — 2.0. 
When Vi is relatively small, the local moment is mainly 
screened by the s-channel conduction electrons, corre- 
sponding to the s-channel Kondo effect. Figure 8(a) 
shows that entropy of the impurity site S'imp is not af- 
fected so much by the p-wave conduction electrons for 
Vi < 0.16. Spin susceptibility Txs shown in Fig. 9 
reaches close to 0.25 of the local moment limit because of 
strong U. At the critical Vi{= 0.1742), the energy spec- 
tra are characterized by the 2-channel Kondo fixed point 
and accordingly the plateau of entropy at (l/2)A:Blog2 
appears at low temperatures. When Vi is larger than the 
critical value, a different type of Kondo effect takes place 
through screening by the p-wave conduction electrons as 
shown in Fig. 8(b). With further increase of Vi, effective 
Coulomb interaction at the impurity site is suppressed 
by the electron-phonon coupling. This effect of renor- 
malization of effective Coulomb interaction appears as 
the higher Kondo temperatures in Fig. 8(b) and also as 
the suppression of Txs in Fig. 9. 

At the 2-channel Kondo fixed point, the thermal aver- 
age of square of ion displacement (g^) converges to 0.739 
at low temperatures. When we substitute the number 
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Fig. 10. (Color online) Temperature dependence of entropy for 
various phonon frequency uj with Vo = Vi = 0.2 and U = 0. 
When Ld is small, Simp shows a plateau at kg log 2. For uj < 0.28, 
the plateau due to the PD is suppressed in the same temperature 
dependence as the usual Kondo effect. Results on the entropy for 
various uj are assembled in the parameter plane of uj and T in the 
lower panel. The red squares are the points where the entropy is 
larger than l.lfce log 2 and the green circles are for 0.9A:b log 2 < 
'S'imp < I-I'sb log 2 and the blue triangles for Simp < O.lfce log 2. 



into (q^), the ratio of Ji.i(q^) to Jo,o is estimated to be 



Jo,o 



1.12. 



(4.6) 



The estimated value close to 1 seems to be reasonable. 

4.5 Non-interacting case (U = 0) 

For discussions on the non-interacting case we fix the 
set of parameters Vq = Vi = 0.2 and U — 0.0. In this 
subsection we use frequency of phonon w as a parame- 
ter to change the strength of electron-phonon coupling, 
which is an alternative to control Vi with a fixed uj. 

Entropy of the impurity site S'imp is shown in Fig. 10(a) 
for various uj. As uj decreases, the plateau of entropy at 
fee log 2 develops owing to the formation of the PD. Fig- 
ure 10(b) summarizes results on the entropy in uj — T 
plane. Typical YAK behaviors are observed in the pa- 
rameter range from oj = 0.238 to 0.280. The middle of 
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Fig. 11. (Color online) Temperature dependence of spin suscep- 
tibility Txs(= 0.25Txc) for Vq = Vi = 0.2 and (7 = 0. When ui 
is lowered, Txs{= 0.25Txc) is suppressed since strong electron- 
phonon coupling increases the binding energy of the PD. 
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Fig. 12. (Color online) Thermal average of square of ion displace- 
ment, (g'^), for Vo = Vi = 0.2 and (7 = 0. When uj becomes 
larger, {q'^) asymptotically comes close to 0.5 which corresponds 
to the zero point fluctuation. 



the blank area between the blue triangle and green circle 
regions corresponds to the Kondo temperature for the 
YAK effect. 

In Fig. 11 the spin (charge) susceptibility of the im- 
purity, Txs (Txc), is shown for various uj. In general 
the relation Txs = 0.25Txc holds in the non-interacting 
cases. When temperature is lower than the binding en- 
ergy of the PD, i?pD, spin and charge fluctuations at the 
impurity site are suppressed due to the formation of the 
PD which consists of the spin singlets of doubly occu- 
pied bonding or anti-bonding states. Therefore, unlike 
the usual spin Kondo effect temperature dependence of 
Txs (Txc) is governed by EpD rather than the Kondo 
temperature of the YAK effect. 

Figure 12 shows temperature dependence of the ther- 
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Fig. 13. (Color online) Temperature dependence of entropy in the 
weak electron-phonon coupling regime, Vq = 2Vi = uj = 0.2. 
When U is increased, a plateau of entropy at fce log 2 appears 
because of formation of the local moment. At low temperatures 
the local moment is screened mainly by the s-wave conduction 
electrons. 



mal average of square of ion displacement, {q^}- For all 
Lu, l^q^^ converges to a constant at low temperatures and 
does not show any changes in the vicinity of the YAK 
temperature. Once the PD is formed the phonon states 
are almost frozen and given by the two distinct coherent 
states which are symmetrically shifted right and left from 
the center. It is interesting to note that (q^) for a;=0.230, 
0.245 and 0.260 show minima in the temperature range 
from 10° to 10"^ which may correspond to -EpD- These 
behaviors may be interpreted as formation of the effec- 
tive double- well potential for the ion displacement. When 
temperature is higher than i?pD, the dynamics of the ion 
is governed by the original harmonic cage potential. On 
the other hand, below the energy scale of -EpD, the ion 
tunnels between the left and right minima corresponding 
to |L) and |R). 

^.6 Effect of finite U in the weak electron-phonon cou- 
pling regime 

Next we choose the set of parameters Vq = 2Vi = 
u = 0.2, where screening by the s-wave conduction elec- 
trons dominates over the p-wave ones. To investigate 
the effects of Coulomb interaction, we show tempera- 
ture dependence of the physical quantities, the entropy 
of the impurity ion 5i„ip in Fig. 13, the spin susceptibil- 
ity Txs in Fig. 14(a) and the charge susceptibility Txc 
in Fig. 14(b). When U is small, S'imp monotonically de- 
cays to zero. Txs and Txc also go to zero without any 
plateau. Therefore, we conclude that the physical prop- 
erties correspond to the RFC regime in Fig. 4. When 
U increases, the impurity orbital begins to behave as a 
local moment. Actually, the plateau region in Fig. 13 is 
broadened with increasing U . Correspondingly, Txs is 
enhanced and shows a maximum close to 0.25 while Txc 
is suppressed, which are signatures of a local moment. As 
already mentioned, the energy spectra are characterized 
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Fig. 14. (Color online) Temperature dependence of Txs shown in 
the upper panel and Txc in the lower one in the weak electron- 
phonon coupling regime, Vq = 2Vi = oj = 0.2. The relation, 
Txs = 0.25Txc, is valid for U = 0. When U is increased, Txs is 
enhanced while Txc is suppressed. 



by the s-type fixed point for any U. However the physi- 
cal properties cross over from the RFC behaviors to the 
s-chK ones with increasing U. 

4-7 Effect of finite U in the strong electron-phonon cou- 
pling regime 

In this subsection we show results of NRG calculations 
for a strong electron-phonon coupling case, Vq = Vi = 
LU — 0.2. In Figs. 15(a) and 16 we see two distinct tem- 
perature dependence of S'imp in the course of release of 
entropy from the plateau at fee log 2. 

For U < 1.4 shown in Fig. 15(a), the electron-phonon 
coupling is relatively more important than the Coulomb 
interaction. In this region, the Coulomb interaction can 
be treated as a perturbation to the PD. 

For the non-interacting case the plateau of entropy at 
fee log 2 is not suppressed even at low temperatures. Be- 
cause of the strong electron-phonon coupling we expect 
that the YAK temperature is much lower than 10~^°. 

For U > 0, the low temperature energy spectrum is 
of the p-type. In §3.2, we showed that the Coulomb in- 
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Fig. 15. (Color online) (a) Temperature dependence of entropy 
Simp for U < 1.2 in the strong electron-phonon coupling regime, 
Vo = Vi = id = 0.2. Plateaus at log 2 originate from the de- 
generacy of the PD. At (7 = the binding energy of the PD is 
so large that Yu and Anderson type Kondo temperature is lower 
than 10"^". For finite U, the degeneracy of the PD is lifted by 
the Coulomb interaction. (b)Temperature dependence of Simp 
for U = 0.01 and its fit by a Shottoky type entropy with the 
excitation gap A = 3.9 x 10~^. 



teraction lifts the degeneracy of the PD in a similar way 
as a transversal magnetic field lifts degeneracy of a free 
spin. When U is smaller than 1.4, it is expected that 
low-energy physics of the impurity is identified as Yu 
and Anderson type Kondo effect under the "transverse 
field" , which is characterized by p-YAK in the phase di- 
agram, Fig. 4. Therefore, it is useful for us to compare 
the present results with those of conventional spin 1/2 
Kondo model under magnetic field. It is shown for the 
latter model that when the magnetic field is larger than 
Kondo temperature behavior of heat capacity approaches 
to a Schottky form (see Fig. 15 in Ref. 39). The lower 
panel of Fig. 15 shows Simp at U = 0.01 fitted by the 
Schottky form with A = 3.9 x 10~^. Let us compare this 
A with the energy gap between the PD, namely ApD. 
To estimate ApD it is appropriate to use a renormalized 
coupling constant A instead of the bare A(= —Vi/ui). 




Fig. 16. (Color online) Temperature dependence of entropy Simp 
for U > 1.4 in the strong electron-phonon coupling regime, 
Vq = Vi = LJ = 0.2. Plateaus at A:Blog2 come from the de- 
generacy of the local moment for the relatively strong U. The 
low temperature energy spectra are characterized by the p-type. 



We can estimate A from the NRG calculations by us- 

ing the relation (^q^^ = (16A -I- l)/2. In this way ApD 
is estimated to be 2.62 x 10^^ for U — 0.01 by using 
(g^) = 5.75. It is important to note that there is a no- 
table deviation at low temperatures from the Schottky 
form. It indicates that energy spectrum at low tempera- 
tures is not gapped but continuous, which is consistent 
with the fact that even with a finite U{y^ 0) the low tem- 
perature properties are described by a local Fermi liquid 
coupled with the p-wave conduction electrons. 

We comment on the transition from the s-YAK regime 
to the p-YAK regime in the case of strong electron- 
phonon coupling. Let us start from the noninteract- 
ing case, C/ = 0. As we have discussed in §4.5 Kondo 
temperature of the YAK effect Tyak is generally small 
and becomes extremely small as the binding energy 
|i?PD| gets larger, see Fig. 10(b). The Kondo energy, 
fceTYAK, is nothing but the energy gain at the s-type 
fixed points. On the other hand, the Coulomb interac- 
tion is favorable to the p-type fixed points and Apo 
is its energy gain. Therefore, the transition takes place 
when fcelYAK ~ ApD. Since fceTVAK is extremely small 
for strong electron-phonon coupling, the transition takes 
place close vicinity of the U = line. 

For U > 1.5 electrons in the impurity orbital behave as 
a local moment and it is quenched dominantly by the p- 
wave conduction electrons. Actually, in Fig. 16 there are 
plateaus of entropy at log 2 and the entropy is released 
in the same temperature dependence as the usual Kondo 
effect. 

Lastly, let us discuss the intermediate parameter re- 
gion, see Fig. 4. In the region between green and blue 
dashed lines, the entropy decays to zero without any 
plateau. It seems that the Coulomb interaction and the 
electron-phonon coupling are competing and almost bal- 
anced each other, leading to almost monotonic temper- 
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Fig. 17. (Color online) Temperature dependence of Txs (a) and 
Txc (b) in the strong electron-phonon coupling regime, Vq = 
Vi = LJ = 0.2. For the [7 = 0, the relation Txs = 0.25Txc is 
valid. Unlike the weak electron-phonon coupling case, the en- 
hancement of Txs by Coulomb interaction is weak for U < 2.0. 
Txs starts to show a maximum for U > 2.0 and becomes close 
to 0.25 as U is increased. 



ature dependence of Txs as shown in Fig. 17(a) and 
Txc in (b). Of course the energy spectra at low tem- 
peratures are always characterized by the p-type except 
for U = 0.0. However the nature of the physical proper- 
ties changes from the PD with a small energy splitting 
due to the Coulomb interaction for U < 1.4 to the local 
moment regime for U > 1.5. The crossover between the 
two regimes is rapid but continuous. 

4-. 8 The line of 8- channel Kondo fixed points 

Finally we discuss physical properties on the line of 
2-channel Kondo fixed points shown in Fig. 4. On the 
line, energy spectra are always characterized by the 2- 
chK type. For the strong Coulomb interaction regime, 
the effective Hamiltonian is given by eq. (4.4), which is 
equivalent to the conventional 2-chK model. It is well 
known that the 2-chK effect is realized when the two 
types of screening channels are fully balanced. In the 
weak correlation regime, however, the degrees of freedom 




Fig. 18. (Color online) U versus (g^) and (S^) with Vq = u) = 
0.2. The sets of the parameters, Vi and U, are chosen to corre- 
spond to each points on the red line with squares in Fig. 4. 



screened by the conduction electrons are not the spin 
degrees of freedom of the impurity orbital but the ion 
displacements of the PD. 

To observe the crossover behavior, we show thermal 
average of square of the ion displacement (g^) and mag- 
nitude of the local moment {Sj') calculated at low tem- 
peratures. Figure 18 shows (g^) and (Sj) calculated for 
the sets of parameters on the line of 2-chK fixed points 
in Fig. 4. When U increases along the line, (g^) is sup- 
pressed by Coulomb interaction and approaches to 0.5 
which corresponds to the zero point fluctuations, (q^) 
decreases monotonically, showing a concave curve as a 
function of U. On the other hand, the magnitude of the 
local moment {S'j') increases from 0.375 to 0.75. The for- 
mer number corresponds to the free orbital limit, while 
the latter to the local moment limit. The curve of {Sj') is 
concave-downward. These results support the view that 
the degree of freedom to be screened by conduction elec- 
trons smoothly transfers from the ion displacement to the 
spin degrees of freedom of local moment with increasing 
U. 

5. Conclusion and future problems 

In this paper we have discussed Kondo effects when a 
magnetic ion is vibrating in a harmonic potential. The 
central issue is the competition between the Coulomb 
and electron-phonon interactions. To treat the problem, 
a generalized impurity Anderson model is derived by tak- 
ing account of vibrational mode of oscillations. Associ- 
ated with phonon absorption or emission, a new chan- 
nel opens for hybridization of the localized orbital with 
conduction electrons. Suppose the localized orbital is de- 
scribed by an .^-spherical wave, then the new channels are 
hybridization with the £±1 partial waves. In the simplest 
case of the s-wave symmetry of the localized orbital, the 
new channel is the mixing with the p-partial waves. 

To proceed further, we have restricted the ion oscilla- 
tion to one dimension, in which case only one component 
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among the three p-partial waves is relevant. With this re- 
duction, the rotational symmetry of the original system is 
lowered to the inversion symmetry. This model still shows 
rich variety of interesting physics which includes dif- 
ferent limiting cases: conventional s-wave Kondo effect, 
phonon-assisted p-wave Kondo effect and Yu-Anderson 
type Kondo effect. 

The model has been analyzed by using the NRG 
scheme and it has been shown that the low-energy fixed 
points belong to one of the two different classes which are 
characterized by the s-type or p-type fixed points. On the 
boundary we find generically the 2-chK fixed points. Sur- 
prisingly, not only the s-chK but also the YAK regions 
are continuously connected to the local Fermi liquid. It 
is also a new finding that the YAK effect which has been 
discussed for a system with strong electron-phonon cou- 
pling but without any effect of electron-electron correla- 
tion is actually very sensitive to the Coulomb interaction. 
With introduction of small U, the fixed point of the YAK 
effect changes from the s-type to the p-type. 

Concerning the Kondo effects of a magnetic ion vi- 
brating in the sea of conduction electrons, there are a 
number of interesting questions to be addressed in fu- 
ture. Response to magnetic fields may be different de- 
pending on where the system is located in the phase di- 
agram. Especially, effect of magnetic field for the YAK 
effect would be very interesting in the region around the 
boundary of the 2-chK. For the present model interesting 
phase transitions take place close to the region of inter- 
mediate electron-phonon coupling, Vq ~ Vi and Vi ~ iv. 
When we consider a strongly anharmonic potential like 
a double well one, the condition for the occurrence of the 
YAK effect would be drastically modified. In view of the 
fact that there are some experimental indications that 
magnetic ions in some skutterudite compounds migrate 
among off-center positions, this might be an interesting 
question relevant to this class of materials. ^""^^^ 
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Appendix: Expansion of the overlap integrals 
within 0{Q) 

In this appendix we derive hybridization terms ex- 
panded within the first order of ion displacement by us- 
ing the generalized gradient theorem. This theorem 
is applicable to the present model because the rotational 
symmetry SO{3) is assumed in the total system including 
both the ion and the conduction electrons. First, let us 
transform the coordinate system from the conventional 
Cartesian form to the spherical tensor form. The position 
vectors of a conduction electron and the ion displacement 



Q are transformed to 

^ r • ^ 

ri,±i = T-^{rxTiry}, ri,o = Tz, 

Qi,±i = T-^{Qx i^Qy), Qi,o = Qz- 

One can show that the set of differential operators 
{Vi,to} defined by {ri.m} are first rank spherical tensor 
operators. Within 0{Q) the overlap integral is written 
as 

V(k;em) (Q) 

J d^re-^'-^iiJemir) - Q ■ VV^„(r)}. 

Q • V is expressed in the spherical tensor form by 
J2mQi,m'^i-,m bccausc of the SO{'d) symmetry. 

The spherical tensors can be expressed in terms of the 
spherical harmonics \i,m) and the Clebsch-Gordan co- 
efficients (^£i,mi;£2,m2\i,m) by 

{i',m'\Wi,rm\i,m) 



i + 



2i + 



^(£,m;l,mrK',m'){|:-^} 



After rewriting e~''" '' by using the spherical wave expan- 
sion, the integrals expanded within the first order of Q 
are given by 

V{k;em) (Q) 

1 

mi — — 1 

where r^^~^^^^) are defined by 



^^—j{i,m;l,mi\e-l,m + mi). (A-4) 

The radial integrals are represented by the spherical 
Bessel functions and the radial parts of t/Jemir) as 

Re{k) = )J^g j drr^ji{kr)ue{r), (A-5) 
n+i{k) 
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I An 



drr^jt+i{kr)\ £ - ^ (A- 6) 




drr ji_i{kr)\ — 



d_ 

dr 



£+1 



ndr). (A.7) 



Lastly we neglect the k dependence of the radial integrals 
and use the values at the Fermi momentum fcp- After 
introducing the dimensionless ion displacement q{= aQ) 

we obtain eq. (2.7), where Ve and V^^^'™ ^ 
Vi = r^ReikF), 



are defined by 
(A-8) 



a 



(A-9) 
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